bo 



Elliptic Orbits with a Non-Newtonian 
Eccentricity 

F.T. Hioe* and David Kuebel 
Department of Physics, St. John Fisher College, Rochester, NY 14618, 

and 



Department of Physics & Astronomy, University of Rochester, Rochester, NY 14627 

August 2, 2012 



Abstract 



o 
cr 

>' | ' It is shown that the lowest order general relativistic correction pro- 

duces elliptic orbits with a non-Newtonian eccentricity. 
PACS numbers: 04.20. Jb, 02.90.+p 

In a weak gravitational field, the general relativistic effect of a massive object 
such as a star that produces a precessional motion [1] to an otherwise Newtonian 
elliptical orbit of a particle (such as a planet) is well known and well noted for 
its historical significance. The precessional angle is approximately given by 
6Tr(GM/hc) 2 , where M is the mass of the star, h is the angular momentum per 
unit rest mass of the particle, G is the universal gravitation constant, and c is the 
speed of light. Like the precessional angle of a particle in a weak gravitational 
field, for small s = GM/hc, most lowest order general relativistic corrections are 
known to be of the order of s 2 and higher. A general relativistic correction of the 
order s, on the other hand, is uncommon and is a principal result that we shall 
present in this Note. Specifically we shall present elliptic orbits with eccentricity 
given by \/2s; that is, we shall present a general relativistic effect of order s that 
makes circular Newtonian orbits elliptical, and the resulting elliptical orbits are 
non-precessing if terms of order s 2 and higher can be neglected. In contrast, two 
new examples for hyperbolic orbits in which the lowest order general relativistic 
corrections are of the order s 2 are also presented. 

We start with one of the analytic solutions for the orbits in the Schwarzschild 
geometry that we presented in our papers [2-6] . We first introduce the parame- 
ters used in the analysis. The massive spherical object (which we call a star) of 
mass M sits at the origin of the coordinate system. Let the coordinates r and 
4> describe the position of the particle relative to the star. If [a;' 1 ] = (t,r,9,4>), 
then the worldline x M (r), where r is the proper time along the path, of a particle 
moving in the equatorial plane — 7r/2, satisfies the 'combined' energy equation 

[1] 
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where the derivative ■ represents d/dr, a = 2GM/c 2 is the Schwarzschild 

radius, h — r 2 (j> is identified as the angular momentum per unit rest mass of the 
particle, and the constant n — Ej (m c 2 ) is identified to be the total energy per 
unit rest energy of the particle, E being the total energy of the particle in its 
orbit and tuq the rest mass of the particle at r = oo. 

By using the dimensionless distance q = r/a of the particle from the star 
measured in units of the Schwarzschild radius and another dimensionless quan- 
tity U defined by 

eq.(l) reduces to the following simple form 
where 



1 2 
52 = T2- S 

g 3 = — + - s 2 --K 2 s 2 = — -—s 2 + -(l-e 2 ) S \ (4) 
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and where 



and 



2 1 h 2 C 2 (K 2 -l) K 2 -l 

6 " 1+ (GM) 2 + ^ 

The use of the dimensionless distance q led naturally to two dimensionless 
parameters k 2 (or e 2 ) and s 2 for characterizing the orbit. As was pointed in 
our previous work [2-6] , the use of the parameter e 2 makes the correspondence 
to the Newtonian case much easier to see. To demonstrate this, we use eqs.(6) 
and (1) to write 



e 2 




GM \ 2GM ■ 2 . 
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and compare this expression with the Newtonian eccentricity ejv, which, us- 
ing the derivative ■ to represent d/dt, t being the ordinary time, can be expressed 
as 



-n 




(8) 



Comparing the two above expressions, it is seen that e — > ejy in the Newto- 
nian limit implies the approximation r — > t and c — > oo. However, since setting 
c = oo is not consistent with reality, we will proceed by stating that e — > ejy if 
we take the approximation r — > t and 

GmY 2GM - 2 

— » — r <b ■ ( 9 ) 
r 2 (f> ) 

We use the coordinates (e 2 ,s 2 ), where — oo < e 2 < +00, < s < 00 of a 
parameter space for characterizing the two regions which we call Regions I and 
II for different types of orbits [5,6]. Region I is mathematically characterized 
by A < and Region II is characterized by A > where A is the discriminant 
of the cubic equation 




4U 3 -g 2 U-g 3 = (10) 



that is defined by 



A = 27^3 ~ 92 (11) 

and where gi and g 3 are defined by eq.(4). For the case A < 0, the three 
roots of the cubic equation (10) are all real. We call the three roots ei,e2,e 3 
and arrange them so that e\ > ei > In this paper, we are interested only in 
the orbit solution for which A < 0, e\ > e 2 > U > e 3 applicable in Region I. 
The equation for the orbit is [2,3] 

- = I+4e 3 + 4(e 2 - e 3 )sn 2 {j(j),k). (12) 
q 6 

The constant 7 appearing in the argument, and the modulus k, of the Ja- 
cobian elliptic functions [7] are given in terms of the three roots of the cubic 
equation (10) by 



7 = (ei-e 3 ) 1/2 , (13) 

k 2 = 6 -^. (14) 

ei - e 3 



where e\, e 3 are given by 
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ei = 2 I — ) cos 



12/ V 3 / 
62 = Hit) COS (- + -)' 



e3 

and where 



3 3 
i/2 /0 2tt\ 
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cos6» = g 3 ^-3 J . (16) 

A typical orbit given by eq.(12) (not on any one of the three boundaries) in 
Region I is a precessional elliptic- type orbit for e 2 < 1, a parabolic- type orbit 
for e 2 = 1, and a hyperbolic- type orbit for e 2 > 1 [5,6]. 

For the elliptic-type orbits (e 2 < 1), the maximum distance r max (the aphe- 
lion) of the particle from the star and the minimum distance r m j n (the perihelion) 
of the particle from the star, or their corresponding dimensionless forms g max 
(= r mM /a) and g min (= r min /a), are obtained from eq.(12) when 7^ = and 
when -f(j) = K(k) respectively, where K(k) is the complete elliptic integral of 
the first kind [7], and they are given by 

^+4e 3 , (17) 



and 

— = i+4e 2 . (18) 

Qmin «5 

The geometric eccentricity e of the orbit is defined in the range < e < 1 

by 

f'max f*min Qmax QVnin e 2 e 3 /-, «\ 

e = = = — — , (19) 

'"max T" ^min Q'max + 9min 1/6 + e 2 + e 3 

using g max and q m [ n given by eqs.(17) and (18). It has been shown in ref.3 
that in the range < e < 1 that e — > e from below as s — > 0, and that e = e for 
all values of s when s = 1. 

The precessional angle A0 is given by 

= - 2tt. (20) 
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The Newtonian correspondence is approached by making s very small. Sub- 
stituting eq.(4) into eq.(16) and expanding in power series in s, we find 



cos 9 = 1 - 2 • 3 3 eV - 2 2 • 3 3 (1 + 9e 2 )s 6 - 2 • 3 5 • 5(1 + 6e 2 )s 8 + ... (21) 
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To obtain a power series in s for 6, the point e 2 = must be done separately. 
For e 2 > 0, we have 



6 = 2 • 3V3es 2 [l + e~ 2 (l + 9e 2 )s 2 + ...] . (22) 

Expanding e\, e2, 7, k 2 , sn(7</>, A:) in eq.(12), e in eq.(19), and A(f> in eq.(20) 
in the power series in s, we find that the orbit equation (12) can be approximated 
for e 2 > and for small s 2 by 

i = 2s 2 {l - e cos[(l - 5)0]}, (23) 
which, in terms of r, gives the approximate orbit equation 

£ = ^{l-*cos[(l -*)*]}, ( 24 ) 
where e, to the order of s 2 , is given by 

e~ e+ (e^ 1 - e 3 )s 2 , (25) 
and where <5, to the order of s 2 , is given by 

(5~3s 2 . (26) 

<5 is related to the precessional angle A(j> given in cq.(20) by Acj) ~ 2nS ~ 
6tts 2 — 6ir[G M / (he)] 2 and it is independent of e (to the order s 2 ) for < 
e < 00. As an example, a general relativistic elliptic-type orbit with e = 0.8, 
s = 0.0176539 has an exact e — 0.80015. The lowest order general relativistic 
corrections (25) and (26) yield e ~ 0.80023 and 6 ~ 0.000935. The corrections 
to e ~ e and the magnitude of 5 are both of the order s 2 . 

Thus if we can ignore terms of order s 2 and higher, we recover the Newtonian 
orbit equation given by 

l - = ^r( l -z™^)- (27) 

The question that can be posed at this point is whether we can define the 
Newtonian limit by stating that it is the general relativistic result for small s if 
we ignore terms of order s 2 and higher. 

To proceed, we note that the case e = is excluded from the expansion 
given by eq.(22) and it is also clear that e 2 = does not satisfy the condition 
given by eq.(9) and must be treated separately. 

For e 2 = 0, the expansion for 9, instead of eq.(22), is now 

6 = VZ^¥s 3 + ^s 2 + (28) 

and the approximate orbit equation still has the form of eq.(23) or (24) with 
the same S given by eq.(26) to the order of s 2 , but with e, instead of eq.(25), 
now given to the order of s 3 by 
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V2s ( 1 + -s 2 ) . (29) 



Thus we have elliptic orbits that precess with the same angle 5 given by 
eq.(26) but with an eccentricity equal to \/2s + 9\/2s 3 /4 to the order s 3 . If we 
ignore terms of order s 2 and higher, the orbit equation becomes 

- = ^(l-\/2scos0), (30) 
r h z 

which is a (non-precessing) elliptical orbit with a non- Newtonian eccentricity 
that is dependent on the speed of light. To the best of our knowledge, this 
general relativistic elliptical orbit with an eccentricity of the order s is new and 
has never been noted by other authors. As an example, an elliptic orbit with 
eccentricity e = 0.0368035 could come from a general relativistic orbit with e = 
[remembering that e is defined by eq.(7) and not eq.(8)] and s — 0.0259843 for 
which the approximation formula (29) gives e ~ 0.0368032 [the first term alone 
gives V2s = 0.0367474] and eq.(26) gives S ~ 0.0020256. The approximation 
e ~ e = holds if we ignore terms of order s. The orbit for e = from general 
relativity becomes a Newtonian circular orbit if we ignore terms of order s, i.e. 
ignoring the second order terms in s 2 and higher order terms is not sufficient to 
get the Newtonian limit for this case. 

It is clear that the entire region characterized by e 2 < or 

.2 / ■ GmV 2GM - 2 

r +[r<f> < ~^-r<j> (31) 

V r 2 4>J c 

is non-Newtonian in character. This includes all circular orbits that occur 
[6] on the curve s 2 = s' 2 for which k 2 = and s = 0, where s' 2 is given by 



_ l- 9e 2_ %/( l + 3e2) 3 

Sl " 27(1 -e 2 ) 2 



from the "vertex" V at (e 2 ,s 2 ) = (—1/3, 1/12) where the innermost stable 
circular orbit (ISCO) occurs, to the origin O at (e 2 ,s 2 ) = (0,0) where the 
circular orbit has an infinite radius. This curve s[ defines a boundary of Region 
I for which the values of e 2 range between —1/3 and and the values for s 2 
range between 1/12 and 0. All circular orbits precess even though the precession 
angle is not observable [8] , and for small s the precession angle is given by 

. , 2 GM 
Aq> ~ Qns ~ bir^ — , 

c l r c 

where r c ~ h 2 /(GM) is the radius of the circular orbit, and A<f> is non- 
zero unless the radius of the circle is infinite which occurs on s = for zero 
gravitational field. 

The values of e 2 along the s[ curve where the circular orbits occur near s = 
are given by 
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2 ^, o '2 

e ~ — isj . 

For small and for s just above Sj inside Region I given by s 2 = sf + (As) 2 , 
it can be shown in the same manner that to the order As we have elliptic orbits 
similar to eq.(30) given by 

1 GM , /-, „ N , N 

- = -T5-l-V2(Aa)cos0. 33 

r hr 

For the parabolic-type orbit (e 2 = 1), e 3 = —1/12 and the initial distance of 
the particle from the star is given from eq.(17) to be q max = oo and eq.(19) gives 
£ = 1. Thus e = 1 and £ = 1 coincide for all values of s. The orbit equation is 
given exactly by 

^ = 4(e 2 + l)sn 2 ( 7 0,fc), (34) 

and for small s values approximately by 
1 GM 



h2 {1 - cos[(l - £)</>]}, (35) 

where 5 is given by eq.(26) and for which the lowest order general relativistic 
correction to the Newtonian case is of the order s 2 . 

For the hyperbolic-type orbit (e 2 > 1), e 3 is less than —1/12 and cq.(17) is 
not applicable. Instead, a particle approaches the star from infinity along an 
incoming asymptote at an angle to the horizontal axis given by [2,3] 



*> = ^' r '(f il^Vl- (36) 

where 7 and k are defined by eqs.(13) and (14), turns counter-clockwise 
about the star to its right on the horizontal axis, and leaves along an outgoing 
asymptote at an angle ^2 given by 

* 2S ^-t, (37, 
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The minimum dimensionless distance g m ; n of the particle from the star is 
still given by eq.(18) as e 2 is still greater than —1/12. 
In the Newtonian limit for small s, ^1 becomes 

*i * cos" 1 ^) = <j> , (38) 
and the complementary angle ^ = 2ir — ^ 2 — <j) also. If we define 

Q GR = *! + *i (39) 

and 
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6 Newton = 20 o , (40) 

the difference A(f> = 9 Newton — @Gfl, can be taken to be an analog of the 
precession angle given by eq.(20) for a hyperbolic orbit, and for small s and to 
the order of s 2 , it was shown [3] to be given by 



A0: 



6tt - 6cj) + 2(2 + e- 2 )\/e 2 - 1 



(41) 



This result is different from the approximation used by Longuski et al. [9]. 
As discussed in ref.9, an experimental test can be carried out to check this result. 

From eq.(18), the minimum distance r m ; n of the particle from the star is 
given approximately by 



1 



GM 



1) 



1 



{e + lf 



(42) 



Equations (41) and (42) show two examples for which the lowest general 
relativistic corrections to the Newtonian case are of the order s 2 and not s for 
e 2 > 1. 

To summarize the above results, for small values of s, the general relativistic 
correction to the Newtonian elliptic orbit is second order in s for e 2 > but is 
first order in s for e 2 < for which the correction is shown to appear in the 
eccentricity of the elliptic orbit. This division gives a new meaning and signif- 
icance to the parameter e 2 . In particular, there exist non-Newtonian elliptic 
orbits of eccentricity y2s given by eq.(30). 
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